Nuclear spins as quantum memory in semiconductor nanostructures 
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We theoretically consider the possibility of using solid state nuclear spins in a semiconductor 
nanostructure environment as long-lived, high-fidelity quantum memory. In particular, we calculate, 
in the limit of a strong applied magnetic field, the fidelity of P donor nuclear spins as a function 
of time in random bath environments of Si and GaAs, and the lifetime of excited intrinsic spins 
in polarized Si and GaAs environments. In the former situation, the nuclear spin dephases due to 
spectral diffusion induced by the dipolar interaction among nuclei in the bath; in the interest of the 
high fidelity requirements necessary for fault-tolerant quantum computing, we consider the initial 
quantum memory decay caused by a non-Mar kovian bath. We calculate this nuclear spin memory 
time in the context of Hahn and Carr-Purcell-Meiboom-Gill refocused spin echoes using a formally 
exact cluster expansion technique which has previously been successful in dealing with electron 
spin dephasing in a solid state nuclear spin bath. With decoherence dominated by transverse 
dephasing (T2), we find it feasible to maintain high fidelity (losses of less than 10 -6 ) quantum 
memory on nuclear spins for times of the order of 100 fis (GaAs:P) and 1-2 ms (natural Si:P) using 
CPMG pulse sequences of just a few (~ 2 — 4) applied pulses. We also consider the complementary 
situation of a central flipped intrinsic nuclear spin in a bath of completely polarized nuclear spins 
where decoherence is caused by the direct flip-flop of the central spin with spins in the bath. Exact 
numerical calculations that include a sufficiently large neighborhood of surrounding nuclei show 
lifetimes on the order of 1-5 ms for both GaAs and natural Si. Our calculated nuclear spin coherence 
times may have significance for solid state quantum computer architectures using localized electron 
spins in semiconductors where nuclear spins have been proposed for quantum memory storage. 

PACS numbers: 03.67.-a; 76.60.Lz; 03.65.Yz; 76.30.-v; 03.67.Lx 



I. INTRODUCTION 

The motivation for developing a solid state quantum 
computer architecture using localized spins as qubits 
arises primarily from the long presumed quantum coher- 
ence times for spins even in the strongly interacting solid 
state environment. In this respect, nuclear spins are ideal 
since both spin relaxation (i.e., T±) and spin coherence 
(i.e., T2) times are very long for nuclear spins, as com- 
pared with electron spins, due to their weak coupling to 
the environment. The application of a strong magnetic 
field further enhances nuclear coherence by suppressing, 
at least, the leading order relaxation and decoherence 
processes caused by direct hyperfmc coupling between 
nuclear spins and any surrounding electron spins due to 
the large mismatch between electron and nuclear spin 
Zeeman energies. Throughout this paper we assume a 
strong applied magnetic field that defines the z-direction 
of our interacting spin system, and we also neglect, some- 
what uncritically, all effects of any direct hyperfine cou- 
pling between electron spins and nuclear spins assuming 
our system to be entirely a nuclear spin system. The ex- 
istence of localized electron spins in the environment will 
further suppress the nuclear spin coherence, and therefore 
our theoretical values for nuclear spin quantum memory 
lifetimes should be taken as upper bounds. 

Solid state nuclear spins as possible qubits and/or 
quantum memory have earlier been discussed in the 
literature^^^ in various contexts although we are not 
aware of any concrete calculations of spectral diffusion- 



induced nuclear spin decoherence using the cluster ex- 
pansion technique (or any other quantitative techniques 
for that matter) along the line we carry out in this work. 
In particular, the vast NMR literature^ (and the more 
recent literature on the prospective use of nuclear spins 
in the quantum information processing^) typically (and 
uncritically) assume exponential Markovian decay of nu- 
clear spin quantum memory whereas we theoretically 
treat the non-Markovian nonexponential decay which 
dominates the initial decoherence of interest to quantum 
memory applications that require high fidelity. We em- 
phasize in this context that we use the term Ti to signify 
the time scale of the initial decay behavior (defined as 
the time in which the extrapolation of this initial decay 
behavior reaches 1/e). This is in sharp contrast to the 
ordinary definition of Ti where an exponential decay is 
implicitly assumed. 

The theoretical work presented in this paper deals 
entirely with the quantum memory of nuclear spins in 
solid state spin quantum computer architectures. In 
particular, we directly calculate solid state nuclear spin 
decoherence in the nuclear spin bath environment, i.e. 
for a Hamiltonian which contains only interaction be- 
tween nuclear spins. Our work is of direct relevance 
to semiconductor-based quantum computer architectures 
where the idea of using nuclear spins to coherently store 
quantum spin memory has been seriously proposed in the 
literatureti^iii^ As such, although we develop a rather 
general theoretical formalism for nuclear spin decoher- 
ence, we concentrate on GaAs quantum dots and phos- 
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phorus doped silicon, Si:P, for applying our theory since 
much of the solid state quantum computer work in the 
literature has focused on these two architectures.— 

The theoretical question we address is simple and well- 
defined: How long is the (high-fidelity) solid state quan- 
tum memory time for nuclear spins in semiconductor 
structures? We answer this question in some detail for 
our specific model and systems, and discuss, using our 
formalism, how one can enhance the nuclear spin quan- 
tum memory (beyond the simple Hahn spin echo situa- 
tion) by using composite pulse sequences to refocus nu- 
clear spins. Our goal is thus a qualitative and quanti- 
tative understanding of nuclear spin quantum memory 
in GaAs quantum dots and Si:P doped structures tak- 
ing into account the complex interacting environment of 
all the other nuclear spins invariably present in the solid 
state environment. 

In considering nuclear spins as solid state quantum 
memory, we discuss two inequivalent (and perhaps com- 
plementary) situations. The usual situation is a bath 
of unpolarized nuclear spins surrounding the "central" 
donor nuclear spin which serves as our quantum mem- 
ory. In this situation, spin relaxation processes (Ti) are 
suppressed by applying a strong magnetic field assuming 
there is a significant mismatch between magnetic mo- 
ments of the donor and intrinsic nuclei. Instead, the cen- 
tral spin memory dephases (T2) primarily through the 
spectral diffusion mechanism that is induced by the sur- 
rounding nuclear spin bath. In the spectral diffusion 
process for the high-fidelity (initial decay) regime that 
is our focus, the central spin quantum memory is effec- 
tively subjected to a non-Markovian temporally fluctuat- 
ing magnetic field due to the nuclear spin flip-flops in the 
surrounding bath that result from dipolar interactions 
among the bath nuclei. The spectral diffusion, examined 
in Sec. HH sets an ultimate limit to the coherence time for 
nuclear quantum memory storage in an unpolarized bath. 
There are, however, well-established procedures, going 
back fifty years;2ii£ which can enhance (or more accu- 
rately restore) the quantum coherence over longer times. 
These coherence-enhancement protocols involve intricate 
refocusing of the nuclear spin dynamics through well- 
designed composite pulse sequences. We discuss a par- 
ticular kind of coherence enhancement protocol, namely 
the Carr-Purcell-Meiboom-Gill (CPMG) composite pulse 
sequence^ in the context of spectral diffusion induced nu- 
clear spin memory decoherence in Sec. IHII of this paper. 
In Sec. lIVl we consider the complementary situation with 
quantum memory stored in a "central" intrinsic spin in a 
bath that is completely spin-polarized (perhaps through 
a dynamic nuclear polarization technique). This some- 
what artificial situation (because complete nuclear spin 
polarization of the surrounding bath is highly improba- 
ble) is instructive because spectral diffusion is, by design, 
completely suppressed in a fully spin polarized nuclear 
spin bath (having no spin degrees of freedom available to 
cause fluctuations); relaxation, however, is now possible 
via direct dipolar flip-flops between the central nucleus 



and like nuclei (those with the same magnetic moment as 
the central nucleus and are thus still allowed to flip-flop in 
the limit of a strong applied magnetic field) in the bath. 
We solve this problem with exact quantum simulation 
(greatly simplified because of the complete polarization 
of the bath) and discuss the results in Sec. IIV1 We con- 
clude in Sec. [V] with a discussion of the implications of 
our results for solid state quantum information process- 
ing and of the various limitations of our theory. 



II. SPECTRAL DIFFUSION OF A DONOR 
NUCLEUS 

Nuclear storage of quantum information in a solid state 
environment is most naturally placed on donor nuclei 
that are easily distinguishable from the surrounding in- 
trinsic nuclei. It is imperative that the memory is stored 
in a nucleus which is distinct from the surrounding nu- 
clei in some manner so that the stored information can 
be recovered. Several quantum computing architecture 
proposals 1 ^^^ exploit the long-term quantum informa- 
tion storage capabilities which donor nuclei spins can pos- 
sess. In this section, we present theoretical calculations 
of the T2 dephasing of donor nuclear spins in two solid 
state environments of interest for quantum computing: 
Si:P and GaAs:P. Specifically, we present coherence ver- 
sus time information in the context of simple Hahn echo 
refocusing; in the next section, we address multiple pulse 
Carr-Purcell-Meiboom-Gill (CPMG) refocusing. 

Fault-tolerant quantum computation demands high fi- 
delity; quantum error correction typically requires 10 to 
10 6 coherent operations 1 ^ and therefore quantum mem- 
ory loss should be kept below 1CT 4 to 10 -6 . The deco- 
herence of nuclear spins has been well-studied in nuclear 
magnetic resonance (NMR) going back 50 year o 6 ' and 
more recently However, these studies were not moti- 
vated by the need for such high fidelity. In NMR, it is typ- 
ical and appropriate to assume Markovian-type decoher- 
ence, with an exponential decay form. However, we find 
that the initial 10~ 6 — 10~ 4 part of the decay of nuclear 
spin coherence is non-Markovian (nonexponential). For 
our purposes, then, T2 must take on a slightly different 
meaning than that of NMR. We find that the initial spec- 
tral diffusion decay will be of the form exp [— (t/Tz) 71 ], 
defining T2 as a characteristic time scale for the initial 
decay (it is defined as the extrapolated 1/e time of the 
decay although the decay will not generally match this 
extrapolation because it will eventually cross over into 
the Markovian regime with an exponential form). For 
considerations of high-fidelity quantum memory storage, 
the n exponent in the exponential holds as much rele- 
vance as the T2 time itself. 

We have previously reported 1 ^ on a cluster expansion 
technique that is often applicable to spectral diffusion 
problems, giving formally exact solutions for short times. 
In some systems, formally exact solutions are obtained 
out to times that are longer than the decay time of 
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the system; in such cases, the cluster expansion method 
solves the entire problem for practical purposes because 
the tail of the decay is of no consequence (assuming a 
monotonic decrease in coherence with time). Other sys- 
tems, including the donor nucleus spin decay currently 
being investigated, yield convergent cluster expansion re- 
sults for only the initial part of the decay. For this reason, 
and motivated by the strict high-fidelity requirements of 
quantum error correction, we focus on the initial decay 
and our plots show memory loss as a function of time 
rather than exhibiting full, formally exact, decay curves. 

In our model, we take the limit of a strong applied 
magnetic field which completely suppresses any interac- 
tions that do not conserve Zeeman energy. For simplic- 
ity, we also disregard interactions with any electrons; our 
Hamiltonian is solely composed of interactions between 
nuclei. We consider only the spectral diffusion decoher- 
ence mechanism due to the dipolar flip-flops in the nu- 
clear bath surrounding the central nuclear spin. The bath 
Hamiltonian consists of dipolar flip-flop interactions be- 
tween like nuclei (as well as the diagonal I nz I mz interac- 
tion): 



_1 1 - 3 cos 2 i 



In+Irn- if 7n = lm 

otherwise 



^Inzln 



,(1) 

(2) 



The free evolution Hamiltonian is given by TIa +Hb- 
We represent the Hahn echo sequence by r — ► tt — ► r 
which we take to mean: free evolution for a time r fol- 
lowed by an NMR 7r-pulse rotation of the P donor nuclear 
spin about an axis perpendicular to the magnetic field, 
and then free evolution for another r in time. The to- 
tal time for this pulse sequence is t = It. In order to 
measure the decoherence of the P donor nuclear spin, we 
compute the magnitude of the expectation value of a P 
donor spin that is initialized perpendicular to the mag- 
netic field (for maximal spectral diffusion) and evolves 
according to the Hahn echo sequence. When normalized 
to a maximum value of one, this expectation value, as 
a function of r, defines the Hahn echo envelope. Where 
(...) is the quantum mechanical average over initial bath 
states, this Hahn echo envelope may be expressed as^ 



W„(r) = ([C/ Hahn (r)] t C/+ hn (r)) 



(5) 



where U^ hn (r) are nuclear evolution operators in the 
context of a Hahn echo sequence with an initially up (+) 
or down (— ) electron spin such that 



U?{t)U±(t), 



-%H±t 



"Hb ± - ^ A n I nz , 



(6) 
(7) 

(8) 



Here, I n ± denote the raising and lowering operator for 
the nth nuclear spin and j n its gyromagnetic ratio, R nm 
denotes the distance between nucleus n and m, and 9 nm 
the angle of the vector between them relative to the ap- 
plied magnetic field. Flip-flop interactions between un- 
like spins are suppressed via energy conservation due to 
the applied magnetic field. 

The qubit-bath interactions are also due to dipolar 
coupling. (We refer to the central spin as a qubit to 
distinguish it from the nuclear spin bath.) However, in 
this instance we disregard the flip-flop interactions which 
are suppressed via energy conservation as a result of the 
applied magnetic field and differing gyromagnetic ratios 
between the qubit and spins in the bath. Instead, we 
consider just the I nz I m z term which will contribute to 
dephasing: 



7~Ia ~ ^ A n I nz S z , 

n 

1-3 cos 2 1 



.4 



(3) 
(4) 



where S z is a nuclear spin operator for the P donor nu- 
cleus, 7 D is the gyromagnetic ratio of the donor nucleus, 
R n is the distance of nucleus n for the P donor, and 8 n 
is the angle of the vector from the P donor to nucleus 
n relative to the applied magnetic field. Note that the 
InzImz term of Eq. (fTJ) requires an extra factor of 1/2 to 
account for interchanging n and m that is not required 
in Eq. ©. 



assuming ideal 7r-pulses. Our calculations assume com- 
plete disorder of the initial nuclear spin bath states, ex- 
pected in the large nuclear temperature limit (T ^> mK) 
without dynamic nuclear polarization or some such tech- 
nique to generate order; the quantum mechanical aver- 
age, (...), of Eq. ([5]) then simply averages over all possible 
nuclear states with equal weighting (that is, the trace of 
the nuclear spin operators divided by the number of nu- 
clear spins states). 

Using the formalism that we have laid out in Ref. 0, 
we can expand Eq. (5} in a way that includes contri- 
butions from "clusters" of successively increasing size. 
Defining vc(t) as the solution to v Hahn (T) [Eq. ©] when 
only considering those nuclei in some cluster (set) C, a 
"cluster contribution" is recursively defined and com- 
puted byi^ 



v 'c( T ) = v c(j) 



E 

{Ci} disjoint, 
Ci^0, d C C 



(9) 



subtracting from vc(r) the sum of all products of con- 
tributions from disjoint sets of clusters contained in C. 
Defined in this way, a cluster contribution can only be 
significant when interactions between nuclei in the set 
are significant and no part is isolated from the rest. We 
consider only local dipolar interactions in the current 
work, and thus contributions only arise when the nuclei 
in the set are spatially clustered together; hence we refer 
to these sets as clusters. The fcth order of the cluster 
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expansion for the echo [Eq. (JSJ)] is the: 



,14 



(10) 



{Ci} disjoint, i 
Ci^Q, \C,\<k 



and involves clusters of size no larger than k. To sim- 
plify the calculation, the logarithm of the echo may be 
approximated as 



In 



,(*0 



(r) 



E "feW, 

|C|<fe 



(11) 



and corrections to this approximation, discussed in 
Ref. [LH may be computed to test its accuracy and are 
expected to be 0(1 /N), where N is the number con- 
tributing spins in the bath. 

We have performed such cluster expansion calculations 
to successively approximate Eq. ([5]) for two different sys- 
tems. In both systems, we have a P donor atom with 
7 D = 7p = 1.08 x 10 4 (s G) _1 , and we have chosen the 
applied magnetic field to point along one of the conven- 
tional axes directions (e.g., Sj j [001]). In our figures, we 
plot "memory loss" versus total echo time (2r) where 
we define memory loss as one minus the echo envelope, 
1 — w H ahn( T )! an( i we only show results where the cluster 
expansion is rapidly convergent and the estimated cor- 
rection to Eq. (fTTj) is negligible. We show Hahn echo 
results (as well as CPMG results which will be discussed 
in the next section) for GaAs:P in Fig. [TJ in this system, 
7„ = 4.58,8.16, and 6.42 x 10 3 (s G)" 1 for 75 As, 71 Ga, 
and 69 Ga, respectively, and all differ from 7 P . We addi- 
tionally show Hahn echo results (and CPMG results) for 
Si:P in Figs. Hand [3 7n = 5.31 x 10 3 (s G)- 1 for 29 Si 
which also differs from 7 P . Unlike GaAs, Si has stable 
isotopes ( 28 Si, 30 Si) with zero spin. Among its stable iso- 
topes, only 29 Si, which has a natural abundance of 4.67% 
and a spin of 1/2, has a nonzero spin. Isotopic purifica- 
tion can reduce the amount of 29 Si and thereby diminish 
spectral diffusion caused by the nuclear spin bath. For 
generality, we define / to be the fraction of Si that is the 
29 Si isotope. Figure shows results in a natural Si bath 
(/ = 0.0467), while Fig. Oshows, for comparison, results 
in a bath of Si isotopically purified to / = 0.01. 

As an alternative to the cluster expansion, a lowest 
order expansion of D HJ1B (r) [Eq. ((SJ)] in r reveals that 
w H ahn(7") = 1 — 0(t 4 ), and this simple lowest order ap- 
proximation compares favorably with our cluster expan- 
sion results where Vh^^t) <C 1. We can be slightly more 
sophisticated with this small r approximation by apply- 
ing it within the context of the cluster expansion; noting 
that v' c (t) = 0(r 4 ), Eq. (JUJ) implies that v Haha (r) = 
exp [— 0(t 4 )]. As discussed in previous work,— conver- 
gence for large systems (i.e., a large bath) is attained in 
this way by embedding the r expansion, which is not nec- 
essarily convergent for large systems, within the context 
of the cluster expansion. 

In the lowest order r approximation discussed above, 



the GaAs:P result becomes 
ln[w Ha h„(T)] w 



260 /is 



(12) 



The exact (convergent) results plotted in Fig. Q] do not 
visibly differ from Eq. (|12[) ; therefore this time expansion 
approximation is valid in the region in which the cluster 
expansion converges. By applying this lowest order time 
expansion approximation, as well as the lowest order in 
the cluster expansion itself (involving only pair contribu- 
tions), the Si:P result becomes 



lnkJr)]^-/ 2 (-J-) 
V 1.05 ms/ 



(13) 



The f 2 dependence simply arises from the fact that, in 
this approximation, all contributions are from pairs of 
nuclei. Dotted lines in Figs. [Hand [3] show the lowest or- 
der approximation of Eq. ([T3")) for the Hahn echo (as well 
as the lowest order approximation for the CPMG pulse 
sequence which will be discussed in the next section). 
The exact (convergent) results exhibit a slight disagree- 
ment with the lowest order approximation as the cluster 
expansion nears the point of its divergence. 

In previous studies of electron spin dephasin g 14 i 15 we 
noted good convergence of the cluster expansion for the 
full decay of the echoes as a result of a particularly ap- 
plicable intrabath perturbation that treats the coupling 
between members of the bath as a perturbation in the 
Hamiltonian. In that situation, the hyperfine coupling 
between the electron and the nuclei was much stronger 
than the dipolar coupling among nuclei. In the case of nu- 
clear spin dephasing, however, the coupling to our central 
nucleus is no stronger than the coupling between mem- 
bers of the bath. In fact, the effective smallness of r 
(compared to the inverse energies of the problem) is ap- 
parently the only reason for any cluster expansion conver- 
gence in the current study. As a consequence, the cluster 
expansion convergence is only attainable for short times 
(small values of r) as indicated in Figs. [1][3] However, 
we do obtain convergent and accurate results of the ini- 
tial quantum memory loss in the duration of high fidelity 
qubit memory retention which is of greatest significance 
in the context of quantum computing. 



III. MULTIPLE PULSE REFOCUSING 

It is well known in the NMR community that multi- 
ple pulse sequences may be used to prolong qubit co- 
herence for a much greater time than the simple Hahn 
echo sequence. The simplest of these multiple pulse 
sequences is known as the Carr-Purcell-Meiboom-Gill 
(CPMG) 1 ^ scheme in which the Hahn echo sequence 
is simply repeated: (r — * tt — ► r) n . In our previous 
worki^ on electron spin coherence in semiconductors, in- 
dependently confirmed in Ref. IT?! we reported that low- 
order symmetry-related cancellations occur when an even 
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number of CPMG pulses are applied; this leads to sig- 
nificantly longer coherence times of even-pulse CPMG 
echoes compared with the Hahn echo. Also, increasing 
the number of these pulses will increase the overall co- 
herence time with the assumption that the applied pulses 
are ideal. Consistent with these previous findings, we find 
that donor nucleus coherence is significantly enhanced by 
applying a two-pulse CPMG relative to the single pulse 
Hahn refocusing, and further enhanced, in overall coher- 
ence time, by applying a four-pulse CPMG sequence. We 
now consider the application of multiple pulse echo se- 
quences in prolonging nuclear spin quantum memory in 
semiconductor quantum computer architectures. 

If we apply 2v CPMG pulses, (r -> tt -> r) 2l/ , the total 
echo time for this is t = Avt, and the exact echo envelope 
may be expressed [analogous to Eq. ((SJ] as^ 

where U^ PMG (v, r) are nuclear evolution operators in the 
context of a CPMG echo sequence with an initially up 
(+) or down (— ) electron spin such that 

C(^) = [U±(t)U*(2t)U±(t)] U , (15) 

referring to the definition of given in Eq. 0. The 
angle brackets, (...), of Eq. (|14p represent a quantum me- 
chanical average over initial bath states that we take to 
be in a uniform distribution representing a completely 
disordered bath as we did in Sec. HT] 

The cluster expansion is performed as described in the 
previous section using Eqs. © and (TTTj) with the only 
exception that vc(r) should now represent the solution, 
while considering only the nuclei in cluster C, to Eq. (fT4")l . 
We have performed these cluster expansion calculations 
to successively approximate Eq. (fT4")) for the GaAs:P and 
Si:P systems. We again plot "memory loss" versus to- 
tal echo time (4i/r) where memory loss is defined as 
1 — v cpmg (t). Results for 2 and 4 pulses {y = 1 and 2, 
respectively) in the GaAs:P system are shown in Fig. [1] 
where they are compared with Hahn echo results. Simi- 
larly, results for the natural and isotopically purified Si:P 
systems are, respectively, shown in Figs. [5] and [3] where 
they are also compared with corresponding Hahn echo 
results. A spin coherence enhancement when applying 
two pulses versus one pulse (Hahn echo) is apparent in 
the plots for each of these systems. 

As shown in Ref. [Tsl if we apply a lowest order time 
expansion and an intrabath perturbation (with coupling 
between nuclei as a perturbation in the Hamiltonian) 
approximation in the cluster contribution terms that 
go into the cluster expansion calculation, we find that 
v CPM g(t) ~ exp [-0{v 2 t & )] = exp [-0(z/~ 4 i 6 )]. In this 
lowest order approximation, the general GaAs:P result 
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FIG. 1: Numerical results of nuclear spin quantum memory 
loss for a 31 P donor nucleus that replaces an As atom in bulk 
GaAs. We define memory loss as one minus the normalized 
echo and plot this in a log-log scale as a function of the total 
echo time. The dashed line gives the Hahn echo results and 
the solid lines give CPMG echo results for two and four pulses. 
At some point for each type of echo sequence, the cluster 
expansion fails to converge. 

becomes 

lnK PMG (r)] « -- 2 (^) 

= ( , (16) 

V 780 /is/ ' v ' 

as compared with In [i; Hahn (r)] w — (i/520 /is) 4 [Eq. (fT2^) 
rewritten in terms of t\. The exact (convergent) results 
plotted in Fig.Q]do not visibly differ from Eq. (TToT) ; there- 
fore this approximation is valid in the region in which the 
cluster expansion converges. This short time approxima- 
tion equation may serve as a useful educated guess (es- 
timate) at times beyond cluster expansion convergence 
although the decay is expected to eventually exhibit ex- 
ponential, Markovian behavior. If we do extrapolate 
Eq. (|16p and define T2 as the time in which the extrapo- 
lated echo reaches 1/e, then we have T2 = j/ - 67 x 780 /is 
for even CPMG echoes. This gives a factor of 6.5 increase 
of nuclear spin coherence times relative to the electron 
spin quantum dot coherence times reported in Ref. Il5l . 

In the range of cluster convergence, where we have con- 
fidence in the accuracy of our results for the model that 
we have used, we observe high fidelity memory reten- 
tion with a low loss of 10~ 6 up to 80 — 120 [is for two 
or four-pulse CPMG sequences. Increasing the number 
of CPMG pulses can improve fidelity versus time the- 
oretically; however, in practice one will be limited by 
the finite time required to perform each 7r-pulse and by 
the accumulation of errors resulting from multiple pulses. 
Because our analysis assumes that each pulse is ideal, 
which is unlikely in practice, these theoretical results for 
Ti should thus be regarded as an upper bound on the 
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Si:P in natural Si 

decoherence of P donor nucleus 



system, we similarly obtain 
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FIG. 2: Numerical results of nuclear spin quantum memory 
loss for a 31 P donor nucleus in bulk Si. We define memory loss 
as one minus the normalized echo and plot this in a log-log 
scale as a function of the total echo time. The dashed line 
gives the Hahn echo results and the solid lines give CPMG 
echo results for two and four pulses. Dotted lines give cor- 
responding results, for comparison, obtained from the lowest 
order expansions provided by Eqs. (|13[) and (|17[) . At some 
point for each type of echo sequence, the cluster expansion 
fails to converge. 
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FIG. 3: Equivalent to Fig.[2]except that results are shown for 
Si purified to 1% 29 Si. Lowest order results given by Eqs. (|13[) 
and (|17[) are shown by the dotted lines. Isotopic purification 
enhances coherence as predicted in these equations. 



1h[v C pmg(t)] w -v 2 f 
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(17) 



as compared with In. [u H8 , hn (r)] w -/ 2 (t/2.1 ms) 4 
[Eq. I|13p rewritten in terms of t]. These lowest order 
results are shown as dotted lines in Figs. [2] and [3l The 
exact (convergent) results initially match this lowest or- 
der approximation, but show a slight disagreement as the 
cluster expansion nears the point of its divergence. 

As with Eqs. fT2")) and (|16p . the above equations may 
serve as a useful educated guess (estimate) at times be- 
yond cluster expansion convergence. Initially, at least, 
Figs. [2] and [3] show that Eqs. (fl"3)) and (fT7|) provide con- 
servative estimates. If we do extrapolate these equa- 
tions and define T2 as the time in which the extrap- 
olated echo reaches 1/e, then we have, for natural Si, 
12 ms for even CPMG echoes. For a small 
number of pulses, v *~ 1, this gives about a factor of 5 
increase of nuclear spin coherence times relative to the 
electron spin quantum dot coherence times reported in 
Ref. [TH; this comparison factor increases as we increase 
the number of pulses because electron spin decay timers 
scales with a smaller power of v (^°- 53 ). In the range of 
cluster convergence, where we have confidence in the ac- 
curacy of our results for the model that we have used, we 
observe high fidelity memory retention with a low loss of 
10~ 6 up to 1—2 ms for two or four-pulse CPMG sequences 
in natural Si and nearly up to 4 ms for 1% purified Si. 

By implementing CPMG pulse sequences with just a 
few even number of pulses, high fidelity (with loss below 
10~ 6 ) qubit retention times are theoretically observed on 
the order of 100 /is for GaAs systems and on the order 
of milliseconds for Si:P systems. We emphasize that al- 
though we are unable to achieve convergence beyond the 
initial decay (for reasons discussed at the end of Sec. |TT| 
which affects the accuracy of our extrapolated estimate 
for T2, itself, we accurately estimate the initial-time co- 
herent memory loss (i.e., the loss of the first 10~ 4 — 10~ 6 
fraction of coherence) which is the most important ingre- 
dient for quantum computation considerations. 



IV. INTRINSIC SPIN RELAXATION 



nuclear memory time. 

Applying the lowest order approximation to the Si:P 



In addition to storing quantum information on indi- 
vidual nuclear spins, there is also a proposal^ to store 
quantum information in coherent states of nuclear spin 
ensembles. The proposal suggests transferring the spin 
state of a quantum dot electron onto coherent spin states 
of the nuclei within the dot via resonance. Before the in- 
formation is transferred, the nuclei must be polarized or 
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cooled into a so-called dark state that will only interact 
with electrons of one type of polarization and not the 
other; this can be accomplished by injecting polarized 
electrons into the dot»i& Then through the hyperfine in- 
teractions between the electron and the nuclei, one may 
induce Rabi oscillations via resonant tuning of the exter- 
nal magnetic field that will pass the quantum informa- 
tion back and forth between the electron spin and coher- 
ent states of the nuclear ensemble. In this way, one can 
transfer the quantum information from the electron spin 
to the nuclear spins. The electron can then be moved off 
of the dot, and the nuclear ensemble will store the quan- 
tum state for milliseconds 1 - until it is retrieved back onto 
an electron spin state by the same process using Rabi os- 
cillations. 

The coherent state of the nuclear spin ensemble in such 
a proposal can be destroyed by nuclear spin dynamics 
induced, for example, by dipolar interactions between 
nuclear spins. If an electron is kept on the dot during 
the storage period, differences of hyperfine energies will 
suppress dipolar interactions that lead to decoherence; 
however, these differences of hyperfine energies in them- 
selves result in an even more rapid loss of coherence due 
to developing phase differences of the nuclear spins 
Thus it is best to move the electron off of the dot; then 
the quantum memory time for this proposal will be de- 
termined solely by interactions (e.g., dipolar) among the 
nuclei. Assuming that the quantum state is imprinted 
on a nuclear spin ensemble that is fully polarized ini- 
tially, the stored nuclear state will be a superpositon of 
the fully polarized state and each possible state with only 
one spin oppositely polarized; then the coherent state is 
lost, in particular, as the oppositely polarized nucleus, 
of the various superposition states, is transported off of 
the dot by dipolar flip-flop processes. This situation is 
analyzed for specific quantum dot geometries in Ref. fl9L 

We will study the problem of the decoherence of such a 
nuclear spin ensemble state in a more general and simple 
way that provides a lower bound limit of coherence time 
(in the ideal case) . Instead of considering a specific dot 
geometry and studying decoherence as a process in which 
spins are transported off of the dot, we simply consider 
the spin relaxation of a single nuclear spin in a bath of 
fully polarized spins. This provides a lower bound, or 
pessimistic estimate, of the coherence time of an ensem- 
ble nuclear spin state such as the previous proposal and 
gives a general consideration for similar proposals that 
may arise. This simple analysis also gives a nice theoret- 
ical counterpart to spectral diffusion which is not possible 
in this scenario; in the limit of a strong applied magnetic 
field flip-flop processes cannot occur among spins of the 
polarized bath. 

With this motivation, we study a well-defined problem 
to compute the lifetime of quantum information stored 
in a central "intrinsic" (as opposed to a donor) nuclear 
spin, surrounded by a bath of fully polarized nuclei of the 
same type, and only considering nuclear spin degrees of 
freedom in the limit of a strong applied magnetic field. 



The precise scenario being considered here in the context 
of nuclear spin quantum memory lifetime is the following: 
The central spin in a polarized nuclear spin bath has its 
spin opposite to all the other spins in the bath initially 
and we are interested in finding out how this local spin 
excitation moves (i.e., "diffuses," although the process 
is technically not a diffusive process) away from the cen- 
tral spin due to dipolar coupling with surrounding nuclei. 
The central spin serves as the quantum memory which 
"decays" in time as the state evolves into superpositions 
of excitation states involving spins in the bath. The ef- 
fective Hamiltonian for the dipolar interaction, under the 
influence of a strong magnetic field which effectively con- 
serves polarization via energy conservation, is given by 



(18) 



where the I n operators operate on the nuclear spin at site 
n in the basis in which z is in the direction of the applied 
magnetic field. The dipolar coupling is given by 



_ 1 2f-2^ 3 cos nm 

Onm - ~jll n H3 • 



(19) 



where 6 nm is the angle formed between the applied mag- 
netic field and the bond vector linking the two spins. 
In a uniform applied magnetic field, which we assume, 
the total Zeeman energy will be constant when polariza- 
tion is conserved. Therefore the magnetic field does not 
affect the dynamics beyond determining the preferred z 
direction and justifying the polarization conservation ap- 
proximation of Eq. (|18[) . 

In our initial state, we will assume a fully polarized 
nuclear spin bath. We denote the fully polarized state, 
including the polarization of the central spin, as \l N ). 
We denote the state with the spin of nucleus n low- 
ered once (via the lowering spin operator) from its fully 
polarized state, but the remaining nuclei polarized, as 
\n) = /„_ \^ N ). We can store a qubit in the superposition 
of states \l N ) and |0). The state \l N ) is an eigenstate of 
Eq. (jTSJ) and is therefore stable. The |0) state, however, 
can decay as a result of a flip-flop with some neighboring 
nucleus; in this section, we will determine the lifetime of 
the |0) state for two different types of materials, GaAs 
and Si, and characterize the decay curve. We can solve 
this problem exactly by using a Hilbert space composed 
of the \n) states as a basis (since our Hamiltonian does 
not include any interactions that can take us out of this 
state space) . Conveniently, this state space grows linearly 
with the number of nuclei (rather than the exponential 
growth of the Hilbert space that occurs in general) which 
allows this problem to be solved exactly even when there 
are many neighbors that influence the decay. 

We wish to directly compute the probability that the 
system remains in state |0) after allowing the system to 
freely evolve for a time t: 



P (t) = ||(0|exp (-mt)\0)\\ 



(20) 
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To solve this problem directly, we first compute the ma- 
trix elements of TL in the \n) basis and then we can ex- 
ponentiate this matrix via direct diagonalization. The 
off-diagonal elements of TL are 

(n| H \m) = 2Ib nm , V n £ to, (21) 

where I is the nuclear spin magnitude, assumed to be the 
same for all involved nuclei. To obtain this we have used 

I±\I,m) = m{m±l)\I,m±l) , (22) 

where only I + [I, I — 1) and I_ \I,I) are relevant. We 
also have 

(n\ H |n) = -2 b krnI{I - 6 nk - 6 nm ) (23) 

= -2^2 hml2 KmL ( 24 ) 

The first term in Eq. (|24|) is a constant contribution to 
the diagonal and therefore a constant, commuting con- 
tribution to the Hamiltonian that will not affect proba- 
bilities. We can therefore use an effective Hamiltonian, 
TL' , in which 

(n\H'\m) = 2Ib mn ,Vn^m, (25) 
(n\H'\n) = 4lJ2 b nm- (26) 

We may now solve the problem by successively includ- 
ing increasingly distant neighbors of the central spin into 
the simulated system until convergence is achieved. In a 
Bravais lattice with only one type of nucleus, (n\ TL' \n) 
[Eq. |26|) ] is, by symmetry, a constant for all n and is 
therefore irrelevant (another constant energy in the sys- 
tem that plays no role in the dynamics). In lattices with 
mixed nuclei, this may not be completely irrelevant; how- 
ever, we have performed calculations on real lattice sys- 
tems with and without this part of the Hamiltonian find- 
ing that it has little effect on the results except that it 
artificially slows down the rate of convergence as we in- 
clude increasingly distant neighbors (because it does lead 
to a significant finite size effect). 

Calculated results are shown for pure 29 Si and natural 
Si in Figs. |4] and [5j respectively. The 29 Si has a spin 
of 1/2 and 77 = 5.31 x 10 3 (s G) -1 . The other stable 
isotopes of Si, 28 Si and 30 Si, have no net nuclear spin 
and therefore do not contribute to the decay. Si has 
a diamond lattice structure with a lattice constant of 
5.43 A. In natural Si, there is a 4.67% abundance of 
29 Si and the random placement of these nuclei leads to 
uncertainty in the results; the standard deviation is thus 
shown by the error bars in Fig. [5[ 

Calculated results are shown for GaAs in Fig. [Sj There 
are three distinct problems with the central spin being 
75 As, 69 Ga, or 71 Ga which have respective gyromagnetic 
ratios of 7/ = 4.58, 6.42, and 8.16 x 10 3 (s G)" 1 and 
all have spins of 3/2. The GaAs is in a zinc-blende lat- 
tice structure with a lattice constant of 5.65 A. The 
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FIG. 4: Relaxation of a central Si nuclear spin in a bath 
of polarized isotopically pure 29 Si nuclei in silicon's diamond 
lattice structure in the limit of a strong applied magnetic field 
along two directions with extremal results: B along [001] (or, 
equivalently, along any of the conventional axes, x, y, or 2) 
and B along [111]. The dashed curves give the corresponding 
approximate solutions from Eq. pip . 
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FIG. 5: Relaxation of a central Si nuclear spin in a bath of 
polarized natural Si nuclei in silicon's diamond lattice struc- 
ture in the limit of a strong applied magnetic field along one 
of the conventional axes (e.g., [100]). Error bars indicate the 
standard deviation as a result of random isotopic configura- 
tions (4.67% 29 Si and the remaining 28 Si and 30 Si is spinless 
giving no contribution to the relaxation). The dashed curve 
gives the corresponding approximate solution from Eq. (|31[) . 



75 As occupies one fee lattice while the 69 Ga and 71 Ga 
are randomly distributed in 60.4% and 39.6% respective 
concentrations on the other fee lattice. The uncertainty 
in the random isotopic configuration of the Ga sublattice 
is the cause for standard deviations shown by the error 
bars in Fig. [6l 

To approximate the short time behavior, let us first 
expand to the lowest order in powers of t. First note 
that (0|i?'|0) = so there is not a first power of t in 
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Excited Nuclear Spin Decay 
Bulk GaAs, B II [001] 
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FIG. 6: Relaxation of a central 5 As, Ga, or Ga nuclear 
spin in a bath of polarized GaAs nuclei in GaAs's zinc-blende 
lattice structure in the limit of a strong applied magnetic field 
along one of the conventional axes (e.g., [100]). In this limit, 
the central spin only interacts with like nuclei. The 75 As are 
arranged with certainty on one of the fee lattices, while the 
69 Ga and 71 Ga randomly occupy 60.4% and 39.6% of the other 
fee lattice, respectively. Error bars for the two Ga results in- 
dicate the standard deviation as a result of random isotopic 
configurations. The dashed curves give the corresponding ap- 
proximate solutions from Eq. (|31[) . 

the expansion. In order to obtain the next order (first 
nontrivial order), we have 

<0|(^') 2 |0) = £<0|«»<n|«'|0) (27) 

n 

= 4/ 2 ^6 2 , (28) 

so then 

(0| exp (-i«t)|0) « 1 - 2I 2 b 2 n0 t 2 - (29) 

Squaring this and keeping only terms up to i 2 , we finally 
obtain 

P «l-(2/) 2 ^6y 2 . (30) 

n#0 

To better match the shape of the decay curves, we use 
an exponential in time that is consistent with the above 
approximation: 

Po^exp l-/(2/) 2 ]T& 2 £ 2 . (31) 

In order to effectively average over isotopic configura- 
tions, we have now explicitly inserted the dependence on 
the fraction, /, of nuclei in the lattice (or sublattice) that 
is of the same type as the central spin. This equation is 



indeed a good approximation to the initial decay and is 
exhibited by the dashed curves in Figs. [H O and A 
good characterization of the lifetime is the time in which 
this approximation reaches 1/e which is given by 



f = l/ f(2iyJ2t>lo- (32) 

This also indicates how the lifetime scales with isotopic 
purification (e.g., in Si). Using this to to define the life- 
time, natural Si has a lifetime of 1.7 ms, and for GaAs, 
75 As, 69 Ga, and 71 Ga have lifetimes of 1.1, 0.7, and 
0.6 ms, respectively. Note that while there is a greater 
abundance of 75 As than the other nuclear types in GaAs, 
it has a longer lifetime due to its relatively small gyro- 
magnetic ratio. 

V. CONCLUSION 

We have analyzed in this work two complementary 
theoretical issues regarding the prospect for using nu- 
clear spins in semiconductor materials, specifically GaAs 
and Si (with and without P doping), as long-lived quan- 
tum memory. We find in general that solid state nuclear 
quantum memory time could be one to two orders of 
magnitude longer than the corresponding electron spin 
coherence times. Thus using nuclear spins as quantum 
memory in semiconductor nanostructures may be useful 
in the context of solid state quantum information pro- 
cessing architectures. 

We address two complementary situations of solid state 
nuclear spin coherence, both taken in the limit of a strong 
applied magnetic field: (1) quantum memory stored in a 
donor nuclear spin in a surrounding unpolarized nuclear 
spin bath where the main memory loss (i.e., decoherence) 
occurs through the spectral diffusion mechanism associ- 
ated with the non-Markovian temporally fluctuating ran- 
dom magnetic field created by the dipolar flip-flops in 
the nuclear spin bath (Sees. HT1 and llllj) ; and (2) quantum 
memory stored in an intrinsic spin surrounded by a spin- 
polarized nuclear bath where the spin relaxes through di- 
rect flip-flop interactions with neighboring nuclei of the 
same type (Sec. H*V]) . These are mutually exclusive situ- 
ations; spectral diffusion disregards spin relaxation (jus- 
tified in a strong magnetic field), and a polarized bath 
cannot induce spectral diffusion. 

In case (1), we note that our T 2 dephasing times define 
characteristic time scales for the initial part of the decay 
that is governed by non-Markovian dynamics in the bath. 
Our T2 values of 520 /zs for P in GaAs and 9.7 ms for P in 
natural Si are consistent with the long-time exponential 
decay behavior studied in NMR, 12 but it is important to 
note that we expect a different short time echo decay be- 
havior of the form exp [-(t/T 2 ) 4 ] = exp [-(2t/T 2 ) 4 ] for 
the Hahn sequence. This has consequences for the na- 
ture of high-fidelity quantum memory retention that is 
most important for quantum computing considerations. 



10 



Longer memory times may be achieved in case (1) by 
using composite CPMG pulse sequences; with an even 
number of 2v pulses and total sequence time t = 4vt, the 
echo decay behavior is of the form exp [— ^ 2 (t/to) 6 ] = 
exp [-v- A {t/t Q f] = exp [-i 2 T 4 /(16r^)]. The latter ex- 
pression is relevant to a train of pulses with fixed t; 
because we consider the non-Markovian, nonexponential 
regime of initial decay, there is no single T2 for a train 
of pulses that is independent of r. Our figures show con- 
siderable coherence enhancement with just two and four 
applied pulses. Further enhancement can be achieved by 
applying more pulses, but this has the cost of requir- 
ing more frequent pulses that must be applied precisely, 
which has physical limitations. Isotopic purification can 
further extend coherence in Si by eliminating 29 Si nuclei 
( 28 Si and 30 Si nuclei have no free spin moments); with / 
as the fraction of 29 Si in the bath, we find that T2 oc 1 / V7 
(or Ti oc J -1 / 3 when pair contributions dominate the de- 
cay of the CPMG sequence). In principle, it is possible to 
extend quantum memory time in Si indefinitely by elim- 
inating all 29 Si, but this is not possible in GaAs where 
all isotopes have nonzero spin. 

When quantum memory is stored on a nucleus that 
is identical to the nuclei in the surrounding bath (rather 
than a donor) , relaxation is not suppressed by applying a 
strong magnetic field and spectral diffusion is somewhat 
ambiguous because the concept of the central spin is ill- 



defined. In case (2), we have chosen to study relaxation 
of an intrinsic nucleus in a bath of polarized nuclear spins 
in which spectral diffusion is suppressed. We find that 
this relaxation decay takes the form of exp [— (i/Ti) 2 ] ; we 
find relaxation times in GaAs of Ti = 0.6, 0.7, and 1.1 ms 
for a central spin of 71 Ga, 69 Ga, or 75 As, respectively, and 
relaxation times in Si of T\ = 1.7 ms/^/f with / as the 
fraction of 29 Si (Ti ~ 8 ms for natural Si in particular). 
These results are relevant, for example, to the proposal 5 
of quantum memory storage in nuclear spin ensembles 
within quantum dots. We do not consider loss of mem- 
ory due to moving an electron on and off the quantum 
dot. We also do not treat the precise dot-dependent deco- 
herence in such ensembles that results from spins "leak- 
ing" off the quantum dot; such a consideration is made 
in Ref. Il9t The virtue of our analysis is in its simplicity 
and generality. By considering the rate at which a spin 
moves away from a central nucleus, we provide a lower 
bound for the coherence times in any such nuclear spin 
memory proposal in the ideal case where dipolar inter- 
actions among the nuclei dominate decoherence and the 
bath is polarized (otherwise, spectral diffusion also plays 
a role). We study this well-defined problem using exact 
quantum mechanical calculations (i.e., diagonalizing the 
Hamiltonian). 

This work was supported by ARO-DTO, ARO-LPS, 
and NSA-LPS. 
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